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I 1. Introduction. 

OO 

I The time dependent Hartree-Fock (TDHF) describes, when the number N of particules tends to infinity, the 
time evolution of the density operator of a quantum system of N interacting particules in the mean field 

I approximation, in other words, when the interaction between two particules is of order 1/iV. For this limiting 

Ph' process when N +00, see, for eaxample, the recent works of [ANl], [AN2], [BGGM] , [ES], [FKS],.... 



A solution to this equation (TDHF) is a density operator Ph{t), that is to say a trace class operator in 
H = L^(]R"') (if each of the particules in interaction is moving in IR"), selfadjoint, > 0, of trace equal to 1, 
evolving as a function of t, and depending on a semiclassical parameter h tending to 0. 



The aim of this work is the analysis of such a solution when the parameter h tends to 0. The analysis 
^ ' will, in particular, make precise the relationship to the Vlasov equation. To establish this relationship, we 

. may associate to the operator ph{t) two notions of symbols: the Weyl symbol (studied in this paper), and 

the Wick symbol (studied in a second paper). The Weyl symbol may only be used under rather strong 
. hypotheses on the operator Ph(0), hypotheses that will be weakened in the second part. In the two articles, 

' the convergence of different symbols when h ^ will be a convergence in L^(IR^"). 

Since (TDHF) appears as a limiting process when the number of particules N de particules tend to infinity, 
a natural question is the one of the interchange of the two limits, the one where N tends to infinity, and 
the one where the semiclassical parameter h tends vers 0. It is the subject of the article [PP] of F. Pezzotti 
and M. Pulvirenti, where it is shown that the Weyl symbol of the marginal density operator associated to 
' a particule in a system of N particules, admits an asymptotic expansion in powers of /i, and that the Weyl 

symbol tends towards the symbol of a solution of (TDHF), and that the coefficient of in the asymptotic 
expansion of the symbol admits a limit when TV tends to infinity , and that for j = 0, this limit is a solution 
of the Vlasov equation. We observe that in [PP] les limits are in the sense of 5'(IR^"), while in this work 
and in the second part they are in the sense of L^(IR^"'). See also [P] and [GMP]. 

The two equations (TDHF and Vlasov) describe the time evolution of particle density, the first one in 
quantum mechanics, and the second one in classical mechanics. Both are dependent on two potentials V et 
W. The first one, V, is the external potential to which all particles are submitted, while W is the mutual 
interaction potential. In this work, these two functions V and W will be real valued C°° functions on IR", 
assumed to be bounded as well as all their derivatives. Let us state precisely now these two equations. 

In classical mechanics, the density at time Ms a function v{x,£,,t) > in i^(IR^"), whose integral equals 1. 
The classical average potential at the point x associated to the density v{.,t) is 

(1.1) V,iix,v{.,t))=V{x)+ [ W{x-y)v{y,v,t)dyd7j. 
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The function v verifies the Vlasov equation: 



. . dv,„^^dv ^ dVci{x,v{.,t)) dv _^ 

If, at an initial time, the data u(.,0) is in L^(]R^"'), and if it is > 0, these two properties remain true for all 
t e H, and the integral IR^" de v{.,t) is conserved, (see Braun-Hepp [BH]). 

In quantum mechanics, the density at time t is described by a family of self adjoint trace class operators 
Ph{t) > (/i > 0) on "H = L^(]R"), with trace equals to 1. The average quantum potential at the point x, 
associated to the density operator ph{t), is: 

(1.3) Vg{x, ph{t)) = V{x) + Tr{W,ph{t)), 

where Wx is the multiplication operator by the function y ^ W{x — y). We will denote by Vq{ph{t)) the 
multiplication operator by the function x Vq{x, ph{t)). The density operator ph{t) is trace class, self 
adjoint, > 0, and verifies the (TDHF) equation : 

(1-4) ih-ph{t) = -h^[A,ph{t)] + [VMt)),ph{t)], 

where A is the Laplacian. To make precise the meaning given to the commutators, let us recall the notion of 
classical solution (TDHF) introduced by Bove da Prato Fano ([BdPFl], [BdPF2]). Let us denote by C^{n) 
the space of trace class operators H = L^{JR"'). Denote by V the space of operators A in jC^{'H) such that 
the following limit 

lim 

t^o t 

exists in C^CH). This Umit is denoted by i[A,A]. A cl assical solution of (TDHF) (for a fixed /; > ) is a 
map t ph{t) in C^{Wi,C^{H)f]C{m.,V) verifying (1.4). In [BdPFl] and [BdPF2], it is shown that if phit) 
is a classical solution of (TDHF) and if the operator ph{0) is > 0, then for all t e H, one has ph{t) > for 
all t, and the trace of ph{t) is constant. 

Our first aim is to understand the Weyl symbol of Ph{t), whose definition is recalled in section 2. The Weyl 
calculus establishes a bijection, that depends on a parameter h> 0, between tempered distributions on IR^" 
(called symbols), and the continous operators of iS(IR") in iS'(IR"). For every tempered distribution F on 
]R^", and for every h> 0, let Op^^^\F) be the associated Weyl operator (see section 2). For every continous 
operator A of 5(]R") in 5(]R"), and, for every h > 0, let (7'^^'"\A) be the Weyl symbol of A. Hence, the 
equalities F = (j^^^\A) and A = Op^^^\F) are equivalent. 

This work uses results of C. Rondeaux [R], that studies the Weyl symbols of operators in the Schatten 
classes, and in particular of trace class operators. We will denote by H^"''P(IR^") (1 < < oo) the space of 
functions that belong to L^'(]R^"), as well as all their derivatives of order < m. It is proved in C. Rondeaux 
[R] that every function in 11^^"+^'-'^ (IR^") is the Weyl symbol of a trace class operator. One can find in [R] 
a characterization of trace class operators whose Weyl symbol is in W°°'^{JR^^). It is also shown in [R] 
that a trace class operator A enjoys this propriety if and only if its iterated commutators with the following 
operators 

(1.5) p.(/,) = ^_^ Qj{h)^x, 

are trace class operators in "H = i^(IR"). This characterization is the analogue of the one given by Beals 
[B] for bounded operators with their Weyl symbol in VF°°'°°(IR^"). One may find in [R] examples of trace 
class operators whose Weyl 's symbol is not in Z/^(IR^"), or of functions in L^(IR^") that are not the Weyl 
symbol of a trace class operator. (The first example will be given in a more explicit form in the second part 
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of this work.) These resuhs wiU be recaUed in section 2, by making clear the dependance on the parameter 
h in view of semi-classical applications, more particularly in the proposition 2.1. 



Our first main result is the following. 

Theorem 1.1. Let {ph{t))(h>o) be a family of classical solutions of (TDHF), with V and W in W°°°°{JR"'). 

We assume that, for every h > 0, the operator p/i(0) is self adjoint > 0, and that its trace equals 1, that, for 
every h> 0, its Weyl symbol a^'^^\ph{0)) is in H^°°^(IR^"), and that the family of functions 

(1.6) Fu = (27r/i)-Vr^'(p,(0)) 

is hounded in VF°°^(IR^") independently of h in (0, 1]. Then, for every t e IR, the Weyl symbol of the operator 
ph{t) is in W°°^(IR^"), and the family of functions Uh{-,t) defined by: 

(1.7) Uh{.,t) = {27rh)-^ary\ph{t)) 

is bounded in W°°^{JR^'^) independently of h in (0, 1] and of t in a compact set of JR. 

The hypothesis of positivity of Ph{0) serves to ensure that the trace norm of Ph{t) remains constant. Without 
this hypothesis it is the trace that stays constant, and not the trace norm. This hypothesis is verified if 
Ph{0) = (27r/i)"Op^'^(G), with G > in L^QR?'''), where Op^^ (G) is the anti-Wick operator associated 
to G (see section 2). One finds in Bove da Prato Fano [BdPFl] and [BdPF2] the proof of the fact that, if 
Ph{0) > 0, then ph{t) > for every t. Observe that the positivity hypothesis concerns the operator and not 
the symbol. 

It is therefore natural to give an asymptotic expansion in terms of powers of h, of the function Uh{-,t) of 

(1.7) . The first term will be a solution of the Vlasov equation and the rest will be estimated in the norm. 
One can see in A. Domps, P. L'Eplattenier, P.G. Reinhard et E. Suraud [DLERS] a physical formulation 
of this problem. The successive terms depend on the initial data (1.6). If the data depends on h, without 
possessing an asympotic expansion as powers of h, it is inevitable that the successive terms in the expansion 
de Uh{-,t) depend on h. 

Theorem 1.2. Let {ph(t) be a classical solution of (TDHF) verifying the hypotheses of theorem 1.1. Then 
there exists a sequence of functions {X,t) Uj{X,t,h) on IR^" x H (j >0), depending on the parameter 

h > 0, such that 

- The function t — >■ Uj{.,t,h) is C°° on IR valued in iy°°'^(IR^"). For every multi-index [a, 13), there exists 
a function Ca0{t), bounded on every compact set o/lR, such that: 

(1.8) \\d^dluj{.,t,h)U^^^.^) < Ca0{t) 

for allt €TR and h e (0, 1] . 

- One has, if Fh is the function defined in (1.6): 

(1.9) uo{X,0,h)=Fh{X) ujiX,0,h) = 0. (j>l) 

- The function uo{X,t,h) verifies the Vlasov equation: 

- For every N >1, the function Uh{.,t) defined in (1.1) and the function F^^\.,t,h) defined by: 

N-l 

(1.11) F^^\X,t,h) = J2 h^Uj{X,t,h) 



fe=0 
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verify if h G (0, 1]; 



(1.12) 



\\uh{.,t) - F^^\.,t,h)\\mjR-^) < CN{t)h 



where Cn is a function on H, bounded on the compact sets. 
-For every N > 1, the operator p^j^\t) defined by: 



(1.13) 



p^^\t) = {2nhrOpZy\F^^\.,t,h)) 



(where F^^\.,t, h) is the function of (1.11)), verifies: 



(1.14) 



Ph{t)-p^mc^n)<C{t)h''+\ 



In section 5, we will make precise the construction of the uj {j > 1), and we will prove the theorem. 

In the absence of mutual interaction {W = 0), and if one studies the time evolution of bounded operators 
instead of trace class operators, the time evolution of the Weyl symbol is described by the Egorov theorem 



If one only makes the hypothesis that ph{0) is trace class, the Weyl symbol of pft(O) is well defined, but is 
not necessarily in i^(IR^"), which is the natural class to compare with a solution of the Vlasov equation 
A counterexample is outlined in [R] and will be detailled in the second part of this work. If one wants to 
strongly weaken the hypotheses on the initial data p/i(0), the Weyl calculus is not anymore the right setting 
to establish a link with the Vlasov equation. On the other hand, wc will sec in a second part of this work 
that, under much weaker hypotheses than the ones of theorems 1.1 and 1.2, the Wick symbol of the solution 
can be used for the relationship between the two equations (TDHF) and Vlasov. 



In section 2, wc will recall the results of C. Rondeaux [R] by making precise the dependance on the parameter 
h in view of applications to semi classical analysis. The corresponding results on the composition of symbols 
and the Moyal bracket are stated in section 3. The results of sections 2 et 3 are proved in the appendixes A 
et B. For section 2 and appendix A we use a technique of A. Unterberger [Ul] and [U2]. The sections 4 et 5 
are dedicated to the proofs of theorems 1.1 and 1.2. 

2. Weyl calculus and trace class operators. 

We denote hy'H = L^(IR") and by £^{71) the set of trace class operators in H. This space is a normed space 
with the norm defined by : 



We will denote by M^'"'P(]R^") {1 < p < +oo, m integer > or m = +oo) the space of functions F which 
are in Lp(]R^"), as well as all their derivatives up to order m. 

Weyl calculus sets a bijection between the operators A of 5(]R") in <S'(]R"), thus admitting a distribution 
kernel in 5'(]R^"), and tempered distributions on IR^" (symbols). This bijection may depend on a parameter 
h > 0. For every F in cS'(]R2"), ^e set Opl''^\F) the operator A : 5(]R") 5'(]R") whose distribution 



(see [H], [Rl], [BR]). 



(2.1) 




kernel is: 



(2.2) 
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This relationship, understood in the sense of distributions, may be inverted. We will denote cr™*^^ (A) the 
distribution F (symbol of Wcyl of A) such that A = Op^^^^{F). In view of applications to trace class 
operators we can rewrite (2.2) equivalently when F is in L^(]R^") as: 

(2.3) OpZ^'iF) = {nh)-^ [ F{X)^xhdX, 

where, for X = {x,S,) in IR^", T,xh is the "symmetry" operator defined by: 

(2.4) {^xhf){u) = ex(«-)«/(2ar - u) X = {x,^ G M^". 
If A is trace class, one has 

(2.5) C7™"^'(^)(X) = 2"Tr(^ o llxh) X e IR^n^ 

According to Caldoron-Vaillancourt [CV] (see also Hormander [H]) if F is in W~~(IR2"), then Op^''^\F) 
is bounded in L^(IR"). If F is in L^(]R^") the operator Op^'^^\F) is bounded by (2.3) but not necessarily 
trace class. It is shown in C.Rondeaux [R] that, if F is in VF™'P(IEI^"), (1 < p < oo, m large enough), the 
operator Opl^y\F) is in the Schatten class L^iW). For p = 1, we will precisely state this in proposition 2.1. 
If F is in W"°°'1(IR2")^ one has: 

(2.6) Tr{Opl''y\F) = (27r/i)-" / F{X)dX. 

If F is in VF°°P(IR2n) q ]^ w°°'i{^^'') (p > 1, g > 1, i + i = 1), one has 

(2.7) Tr(Opr^'(F)oOpr^'(G')) = (27r/i)-"y F{X)G{X)dX. 

The left hand side makes sense, since from [R], the two operators under composition are respectively £^{'H) 
and jC^CH), in such a way that their composition is trace class. 



The next proposition mainly due to C. Rondeaux [R], gives sufficient conditions (on the Weyl symbol) for 
an operator to be trace class, and gives a characterization of the set of operators whose Weyl symbol is 
in W°°^(IR^"). This characterization is the analogue of the Beals characterization [B], which concerns the 

symbols in VF°°°°(IR2"). 

We note Pj (h) = j the momentum operators and we note Qj (h) the multiplication by xj . 
For every operator A of <S(1R") in 5'(IR"), we set 

(2.8) {ad P{h))°'{ad Q{h)fA = {ad Pi{h))°'K..{ad Qn{h)f"A. 



Proposition 2.1. a) If F is in i^(lR^") as well as all its derivatives up to order 2n + 2, then, for all h> 0, 
the operator Op^^^\F) is trace class, and: 

(2.9) \\OpZ''{F)\U.^n)<Ch-" E h^\"\+\^\y^dSdf^Fh^^^..y 

\a\ + \l3\<2n+2 

b) If A is a trace class operator, and if, for every multi-index (a,/3) such that \a\ + < 2n + 2, the operator 
{ad P{h))"{ad Q{h))^A is trace class, then the Weyl symbol of A is in L^(IR^") and: 

(2.10) (27r/i)-«||ar^'(A)|Ui(K-) < C ^ /i-(H+l/3|)/2||(ad P(/j))"(ad Q(/j))/3a||£X(„), 

|c<| + |/3|<2n+2 
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where the constant C depends only on n. 

c) the following are equivalent : 

i) A family of operators (A/j)o</(<i is of the form = Op^^^\Fh), where {Fh) is a bounded family of 
functions in W°°'^{m?"). 

a) For every h > 0, the operator A^ is trace class , as well as all iterated commutators of A^ with the 
operators Pj{h) and Qj{h) and , for every {a,p), the family of norms 

(2.11) h^-M-m \\{adPih)nadQ{h)fA\\crin) 
stays bounded when h varies in (0, 1]. 

Parts a) and b) arc proved in C. Rondeaux [R], without the parameter h. Part a) needs only minor 
modifications to introduce this parameter but it seems to us better to give a proof of Part b) in appendix A. 
For the sake of clarity it might be useful to recall the well-known analogue of proposition 2.1 for the bounded 
operators and the symbols in W°°°°{M?"), that is to say, the C alder on- Vaillancourt theorem and the Beals 
characterization. 

Proposition 2.2. a) If F is in L°°(]R^") as well as all derivatives up to order 2n + 2, then, for all h > 0, 
the operator Op^'^^^{F) is bounded in H = L^(IR"), and: 

(2.12) \\opr\F)\\c^n) <c h^^''\+\^\y^\\d:diFUoo^^..y 

|a| + |/3|<2n+2 

b) If A is a bounded operator and if, for all multi-indices {a, 13) such that \a\ + \I3\ < 2n + 2, the operator 
{ad P{h))"{ad Q{h))^A is bounded, then the Weyl symbol of A is in L°°(]R^"), and one has: 

(2.13) ||ar^'(A)|U<»(B,..) <C J2 /i-''"'+""^/'llM P{h)nad Q{h) fA\U^n)- 

\a\ + \l3\<2n+2 



Anti-Wick calculus. The definition of this calculus uses coherent states. By this we mean the family of 
functions 'i'xh in .^^^(IR"), indexed by the parameter X = {x,^) G IR^", and depending on /i > 0, defined by 

(2.14) *x,/,(w) = (7r/i)-"/^e-^^e^«-S-2fe^-« X = (a;,0 e M^". 

These functions, called coherent states, will be used to give examples of operators verifying the hypotheses 
of theorem 1.1, and will be also helpful to prove Proposition 2.1 (appendix A). Fore more properties of the 
coherent states and the anti-Wick calculus, see, for instance [CF], [F], [CRl], [LI], [Rl ], [Ul], [U2]. The two 
fundamental properties are the fact that: 

(2.15) I < ^xh, '^Yh>\= e-i^l^-^l' ll^-x^ll = 1, 
and that, for all / and g mV. 

(2.16) <f,g>= (27r/i)-" / <f,^xh>< ^xh, 9 > dX. 

For every function F in L^(]R^") and for every h > Owe set Op^^ {F) to be the bounded operator in L^{R^) 
such that for all / and g inH: 

(2.17) < Opf^{F)f,g >= {2nh)-^ [ a{X) < f,^xh > < '^xh,9 > dX. 
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If F is in i^(]R^") we see that Op^^{F) is indeed trace class in H, and that: 

(2.18) WOpt'^iFnc^n) < (27r/i)-" / |F(X)|dX 
Moreover one has: 

(2.19) Tr{Ovt^{F)) = (27r/i)-" / F{X)dX. 

If F > 0, the operator Op^^{F) is self adjoint and > 0. The Weyl symbol of the operator Op^^{F) is 
given by : 

(2.20) ar''(Opt'^{F))=ei^F, 

where A is the Laplacian on H^". In fact, the operators T,Yh defined in (2.4) and the operator Pxh of 
orthogonal projection on Vect{^xh) verify: 

(2.21) Tr{Pxh^Yh) = e-^^. 
3. Basic Facts on the Moyal bracket. 

The composition of symbols in the Weyl calculus is a very classical field (see Hormander [H ] or D. Robert 
[Rl] for the dependance on the semiclassical parameter h). We need to adapt this to the classes of C. 
Rondeaux symbols, by making precise the dependance on the parameter h. 

We define a differential operator cr(Vi, V2) on IR^" x H^" by: 

(3.1) C7(Vi,V2) = ^ 



dyMi dxjdr]j 
where {x,^,y,r]) denotes the variable of IR^" x IR^". 

Theorem 3.1. For all functions F in W°°2'(IR2") and G in H^°°9(]R2n)^ (p>l,q>l,^ + \ = l),for all 
h> 0, there exists a function Mh{F, G, .) in W°°^{M'^") (Moyal bracket) such that: 

(3.2) [Opr'\n OpZ^'iG)] = OpZ^\Mh{F, G, .)). 
For all integers N >2, one has: 

N-l 

(3.3) Mh{F, G,X)=Y^ h'^CkiF, G, X) + Ri''\F, G, X), 
where the function Ck{F,G,X) is defined by: 

(3.4) CkiF, G, X) = ^^i^ [a(Vi, V2)'=(i^ G)(X, X) - a(Vi, V2)'=(G ® F){X, X) 



and where the function {F, G,.) is in W°°^ (H^") . We can write also ( 3.3) for N = 1, with the convention 
that the sum vanishes, and that i?^^^ = Mh- For every integer £, there exists a constant C such that: 

(3.5) h"^\\V'R^^\F,G,.)U.fj^.n)<C /i^"+^^/'l|V"F||L.(ni-) ||V^G|U,(n,2n). 



|a| + |/3|<«+4n+2+2iV 
|a|>N,|/3|>N 



The operator R^^\F,G)=Opl''y\R!i\F,G,.)) verifies: 

(3.6) ||i?f)(F,G)|Ui(«)<C/i-" /^("+^'/'||V"f^|U.(n,..) ||V^G|U,(K-). 

|c.| + |S|<6n+4+2N 
|a|>N,|^|>N 

This theorem will be proved in the appendix B. We will also use in the second part of this work the well 
known analogue of theorem 3.1, which we recall here in order to allow for its application when needed. 

Theorem 3.2. With the notations of theorem 3.1, if the functions F et G are in W°°°°{^?'^), the function 
r\^\f,G, .) defined by the equality () verify, for any i > 0; 

(3.7) h'^^V'Ri''\F,G,.)Uoo^^.n^<C Yl /i'^'+'^/'llV^FIUoofR^n) ||V'=G|Uoo(n,.„) 

i>N,k>N 
e+2N<j+k<e+2N+4n+2 

The operator r\^^ {F, G) verify: 

(3.8) \\R[''\F,G)\\cin)<C ^ h^'+'^^^WV^ F^^^^^^^ || V'=G|Uoo(n,.n) 

j>N,k>N 
2Ar<j + fc<2JV+6n+4 

4. Egorov's theorem for trace class operators and proof of theorem 1.1. 

We are going to adapt to the case of symbols in L^(IR^") and trace class operators the idea of the proof of 
Egorov's theorem contained in the book of D. Robert [Rl]. The difference with [Rl] comes from the fact the 
classes of operators considered here are the classes introduced by C. Rondeaux and that Hamiltonians are 
time dependent. 

Wc consider a function (.x,t) V{x,t) on K" x K, real valued, depending in a C°° manner on a;, and 
continuously on t. We suppose that, for every a, there exists Gq- > such that: 

(4.1) \d^V{x,t)\<Cc, (a;,i) elR" xH. 
We set: 

(4.2) H{x,^,t) = \^\'' + V{x,t). 
We denote by V{t) the multiplication by V{.,t). We set: 

(4.3) Hh{t) = -h''A + V{t). 

Therefore, Hh{t) = Opl''y\H{.,t)). Let us recall now some facts on unitary propagators ([RS]). 

Proposition 4.1. For all t e IR, let V{.,t) a C°° function on M", verifying (4-1), and depending in a G^ 
manner on t. Let H}i{t) be the operator defined in (4-3). For every f in <S(]R"), and for every s € IR, there 
exists a function denoted by t ^ Uh{t, s)f that verifies: 

(4.4) ih^Uhit, s)f = {Hh{t))Uh{t, s)f, Uh{s, s)f = f. 

The operator Uh{t,s) maps <S(IR") into itself and, by duality, <S'(IR") into itself. One has Uh{s,t) = 
Uh{t,s)~^. One also has: 

(4.5) ih^Uh{t, s) = -Uh{t, s){Hh{s)). 
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For every operator A o/<S(IR") into <S'(]R"'), let us set: 

(4.6) Gh{t, s){A) = Uh{t, s)oAo Uh{s, t). 

One has 



(4.7) 



ih^Gh{t,s){A) = [Hh{t),Gh{t,s){A) 



Gh{s,s){A)=A. 



The analogue of Egorov's theorem for operators in the class of [R] is the following. 

Theorem 4.2. Let F be a function defined on VF°°1(IR2"). Let Ah = Op'^'''"\A). Then, for every t e M, 
there exists a function Fht in W°°^{Wi^") such that: 

(4.8) Gh{t,0)iAh)=Op';:'^\Fht). 

If the function F and the potential V{.,t) depends on a parameter X, while staying bounded respectively in 
W°°^ (TR^") and in VF°°°°(IR^") independently of X, (and oft for V), then the function Fht remains bounded 
in M^°°^(]R^") independently of X, of h in (0,1], and oft in a compact set o/M. 

Following the idea of [Rl], which is related in some sense to Dyson's series, we will express our solution 
Gh{t,0){Ah,) under the form: 

N-l 

(4.9) Gh{t, 0){Ah) = Yl Opr' t)) + h^E^it, h), 

fe=0 

where the functions -Dfc(., t) are in H^°°^(]R^") and -Bjv(t, h) will be a trace class operator with bounded trace 
norm. In a second step, wc will show that the commutators of EN{t, h) with the position and momentum 
operators are trace class operators, and we will estimate their traces. Finally we will rely on the character- 
ization of C. Rondeaux (recalled in proposition 2.1) to show that Gh{t,0){Ah) is itself a pseudodifferential 
operator, with a symbol in 1^°°^(IR^"). 

The construction of the terms Dk{.,t) will use the Hamiltonian flow of H{.,t). For every function G in 
W°°1(IR2"), we call $ts(G) the function on R^n defined by 

(4.10) = = 

Under hypothcsc (4.1), one knows that if G = G\ where {Gx)\^e is a family of bounded functions in 
VF°°i(]R2n-)^ then <^ts{Gx) stays bounded W°°^QR^'') when {t,s) varies in a compact set of IR^ and A in E. 

Lemme 4.3. For every function G in W°°^{M?'") and for {t,s) in JR^, one has 

(4.11) Gh{t,s)(0P;:'y'{G)) = Opr^'i'^tsiG)) + h J'^ Gh{t,t,)(Opr''{R{;ti,s,h))) dtu 

where the function R{.,ti,s, h) is in W°°^(]R^"). IfG belongs to a family of bounded functions in 1^°°^(IR^"), 
it is also the case for R{.,ti,s,h) when {ti,s) lies in a compact set o/M^ and h is in (0,1]. 

Proof of the lemma. Prom definition (4.10): 

^Opr\^ts{G)) = OpZ'\{H{.,t),^ts{G)}). 
9 



With the notations of theorem 3.1, with N = 2, one may write: 

[H,{t),Opr\<i>,,,{G))] = %r\{H{.,t),^,,{G)})+OPr'[R^^\H{.,t) ,$*.(G))). 
Consequently: 

On the other hand, 

^GH{t,s){OPr'{G))-lmt),GH{t,s){OPr''{G))]=0 

By combining these two equalities, noting that for t = s, the two operators Gh{s, s){Op^^^\G)) and 
OPh^^^^ssiG)) are equal, and then by using Duhamel ' s principle, we obtain (4.11), with: 

R{., tus,h) = - ii), $ti.(G)) . 

It is well known that, when F{x,£,) = one has r'^^ (^F, G j = for every function G. Hence: 

R{., h,s, h) = -^R^h^ (y{., h), <^tAG)) . 

By hypothesis, V{.,ti) is in TF°°^(IR") and is bounded, independently of ti. We have seen that ^usiG) is 
in VK°°^(]R^"), independently bounded of ti and of s when {ti, s) varies in a compact set of IR^. According to 
theorem 3.1, applied to the case = 2, it follows that ti,s, h) is in VF°°^(IR^"), independently bounded 
of t\ and of s when {t\,s) varies in a compact set of and h in (0, 1]. 

Proof of theorem 4.2. First step. Let F be a function in H^~1(IR2"). Let Ah = Opl''y\F). Let ^ts{G) be 
the function verifying (4.10). For every t e H, we define a function £'o(t)(.) in W^°°^(]R^") by 

(4.12) Do{.,t) = ^tfi{F)- 

We have seen that this function is in VK°°^(]R^"), independently bounded of t on every compact set of H. 
By lemma 4.3 applied to G = F and s = 0, and from (4.12), one has: 

Gh{t,Q){Ah) = OpZ^' (Do {t)) + h j'^ Gh{tM) {OpT'' {Ri{-M,h))) dti , 

where Ri{., t\,h) stays bounded in W°°^^^^"''> when ti belongs to a compact set of IR and h is in (0, 1]. We 
reiterate, by applying 4.3 with s = ti and G = Ri{.,t\,h). We obtain: 

Gh{t, t^)(OpZy\R^{., t,,h)) = OpZ^' U{t, t,){R,{., t,,h))) +h f Gh{t, t2){OpZ''{R2{; t2, tl,h)))dt2, 

where R2{.,t2,ti,h) stays bounded in W^°°^(IR^"') when (ii,t2) belongs to a compact set of IR^ and h is in 
(0, 1]. We define a function Di{.,t) in W^°°i(]R2n) ^y. 

Di{.,t)= f ^t,ti){Ri{.,ti,h))dti. 
Jo 

10 



This function is in W°°^{JR^"'), independently bounded of t on very compact set of JR. One has, if t > 0: 

Gh{t,0){Ah)=OpZ'''(Do{t)+hD,{t)] f Gh{tM]iOpZy\R2{.MM,h)))dti dt2. 

J Q<tl<t2<t 

By reiterating this process, we obtain, for every A'', the equahty (4.9), with: 

(4.13) EN{t,h)= f Gh{t,tN)(Opl'y\RN{.,tN,...,tuh))] dti...dtN, 

JAjv(t,0) ^ ^ 

where Ajv(t, s) is the set defined, if s < by: 

(4.14) /^N{t,s)={(h,...tN)&M,^, S<h< ...<tN <t}, 

and in a symmetric way if s > i. In (4.9), the Dj{.,t,h) [j > 0), and RN{.,tN, ■■■,ti,h), are in H^°°^(]R^"), 
independently bounded of h in (0, 1], of (ti, ...tjv) in Ajv(i, 0), and of Hn a compact set of IR. 

The second step will consist in getting upper bounds for trace norms of iterated commutators of EN{t,h) 
with the position and momentum operators, with the help of a less precise technique than the one used in 
the first step for powers of h. However this lesser precision will be compensated by the presence of the factor 
in (4.9). In order to do that, we will use the following lemma, that will be useful also in section 5 and 
in a second article, (part II of this work). If an operator A is bounded in L^(IR"), as well as all its iterated 
commutators up to order m, we set: 

(4.15) ir°°{A)= Yl \\{adQ{h)f{adP{h)rAU^^y 

|a| + |/3|<m 

If an operator A in L^(IR") is trace class, as well as all its iterated commutators up to order m, we set: 

(4.16) C*'^(^)= E \\iadQih)f{adPih)rA\U^^^y 

|a| + |/3|<m 

The aim of the next lemma is to show that these properties are stable par the mapping Gh{t, s). 

Lemma 4.4. Let Hh{t) be the operator defined in (4-3), where V{.,t) verifies (i-l). Let Uh{t,s) denote the 
unitary propagator and Gh{t, s) the map of proposition 4-1- Let A be a trace class operator inH = L^(IR"), 
as well as all iterated commutators {ad Q(h)Y{ad P{h))"A for \a\ + < m. Then, for all s and t in IR, 
the operator Gh{t, s){A) is also trace class, as well as all iterated commutators with the Pj{h) and Qj{h) up 
to order m. Moreover, for every compact set K o/M, there exists Ck > such that: 

(4.17) 7r'*''(G^(t,s)(A))<(7x/r*'(^) {s,t)eK\ he {0,1]. 

An identical result holds for bounded operators and for the norms 

We see why this lemma alone does not allow to prove theorem 4.2. One will need, for that, that in the 
expression (4.16), the terms corresponding to the multi-index (a,/3) be multiplied by 

Proof of Lemma. By (4.7) one checks that for every operator A verifying the hypothesis of the lemma, and 
for each of the momentum operators Pj{h), the following equality: 

^^[Pj{h), G^t, s){A)] - 1 [Hh{t), [Pj{h), Gn{t, s){A)]] = 1 [[P^{h), Hn{t)], Gn{t, s){A)] . 
Then it results, by the Duhamel principle, the following equality: 

[Pj{h),Gh{t,s){A)]=Gh{t,s)([Pj{h),A]) + ^ j'^ Gh{tM)(j[Pj^^^ 
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We have an analogous equality for the position operators Qj{h). One has: 

[Pj{h),H,{t)] = mh),HH{t)] = 2zhP,{h). 

We therefore deduce: 

[Pj{h), Gh{t, s){A)] = Gh{t, s){[P, ih),A]) - ^* Gh{t, ti) ( f^Q^^'\ Gh{ti,s){A)] ) dt 



[Qj{h),Gh{t,s)iA)] = Gh{t,s)(^[Qj(h),A]^ +2 J Gh{t,ti)(^[Pj{h),Gh{ti,s){A)\yti. 

If A and its commutators with Pj{h) and Qj{h) are trace class, we first observe that [Pj{h),Gh{t,s){A}] 
is a trace class operator since Gh(t,s) conserves L^(Ti.). Reporting in the second equality, we sec that 
[Qj{h), Gh{t, s){A)] is itself a trace class operator, and that the upper bound (4.17) is proved for m — 1. We 
pursue the same reasoning to prove (4.17), by induction, for every m. The analogue (4.17) for the bounded 
operators is proved similarly. 

Proof of theorem 4-2. Second step . Following proposition 2.1, it suffices to show that, for every multi-index 
(a, /3) and for every compact set K of IR, there exists Cafix > such that: 

(4.18) h"-^\"\+\^^{ad P{h)r{ad Q{h)fGh{t,0){Ah)\\cHH) < Ca^K, 

at dz K and h G (0, 1]. In order to achieve this, one will use the asymptotic expansion (4.9) up to an order N 
that will depend on a and /3. Since from the first step the Dj{.,t,h) {j > 0) of (4.9) belong to TF°°1(]R2")^ 
and are bounded independently of h £ (0, 1] and of t in a compact set of K, the proposition 2.1 shows that: 

P(h)r{ad Q{h)fOpr'^\D^{.,t,h))\\ciin) < Ca^K 

if t G K et h £ {0, 1]. Let us derive now an upper bound that is the analogue for the term E]\[{t, s, h). For 
this, we use the form (4.13) of E^it, s, h), and we apply lemma 4.4 with {t, s) remplaced by {t, t^) and A by 
Op^^^\RNi-,tN, ■■■,ti,h)). Since i?jv(., iiv, ti, /i) is in W°°^(Si^") and is independently bounded for h in 
(0, 1], of {ti, ...^Ar) in Ajv(0,i), and for t in a compact set of IR, proposition 2.1 shows that for every integer 
m > 0, for every compact set K of IR, there exists C > such that : 

h"ir'''{Opr'\RN{;t^,.:,h,h))) < CmK 

if < /i < 1, (fi, ...tiv) € Aiv(0, t), and t G K. Hence by lemma 4.4, we deduce that the iterated commutators 

of Gh{t,t]\r){Op^'^'^''{RN{-,tN, ■■■,ti,h))) with the position and momentum operators are themselves trace 
class, and that there exists another constant CmK such that: 



h^i;:''''(GH{t,tN){Opr''{RN{.,tN,...,tuh)))) < Grr 

We can therefore write, if Ah = Opl''y\F), for every multi-index (a, /3) and for every integer A'': 

hP\\{ad Q{h)f{ad P{h)rEN{t,h)\\cHU) < C^^n 

By reporting this in (4.9), and by choosing A'' = |a| -|- |/3|, one deduces (4.18). Using the characterization of 
C. Rondeaux (proposition 2.1), the theorem 4.2 follows. 

□ 
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Proof of theorem 1.1. 

Let Ph{t) be a classical solution of (TDHF) verifying the hypotheses of theorem 1.1. Let us denote by Vh{t) 
the operator of multiplication by the following function 

(4.19) x^V^{x,ph{t)) = V{x)+Tr{W,oph{t)), W^{y) = W{x - y). 

Under the hypotheses of theorem 1.1, p,,(0) > 0. By the results of [BdPFl] and [BdPF2], ph{t) > for all 
t, and the trace norm of ph{t) is constant. Since all the derivatives ofVetW are bounded , one has: 

(4.20) \d^Vg{ph{t)){x)\ < C„, {x,t) e B" X IR. 

Let Hi^{t) denote the operator defined in (4.3), where V{t) is the multiplication by Vq{x, ph{t)). With these 
notations, the equation (TDHF) has the form: 

(4.21) ih^E^ = [H^{t),pnit)]. 

We note that V{t) depends on h. However in theorem 4.2, V{t) may depend on a parameter that could be 
h. The only requirement is that Vq{., Ph{t)) be bounded in 1^°°°°(]R") independently of h, which is the case. 
By denoting Gh{t, s) the unitary propagator associated to the Hamiltonian Hh{t) as in proposition 4.1, one 
hcis tliGrGforG* 

p^{t) = GHit,0)ipHm = Gh{t,0){OpZ^\Fn)). 
Theorem 1.1 is therefore a particular case of theorem 4.2. 

5. Proof of theorem 1.2. 

We are going to state precisely the "explicit" construction of an "approximate solution to order A''", denoted 
by p\l^\t), of the (TDHF) equation. The exact solution ph{t) is determined by the interaction potentials 
V ct W, which belong to V[^°°°°(]R"), and the initial data Ph{0)- Under the hypotheses of theorem 1.1, the 
operator ph{0) is > 0, and is of the form 

p,(0) = {2nhr OpZ'\f,) 

with Fh in W°°^{M'^"'). This function may depend on h, but stays bounded in W^^iM"^") independently of 
h in (0,1]. 

We look for an approximate solution with the ansatz (1.13), where F^^\t, h) is a function IR^" of the form 
(1.11), the functions Uj{.,t) being in W°°^(1R^"). We associate to this approximate solution the average 
quantum potential: 

(5.1) Vqix,pi''\t)) = Vix)+TriwJ,f\t)). 
One knows that, if F is in M^°°1(IR2») and G in M^°°°°(IR2"), one has: 

Tr{Op';;^''\F) o Opl''y\G)) = (27r/i)-" / F{X)G{X)dX. 

Therefore, if p')^\t) is defined by (1.13), we have: 

(5.2) Vq{x,p^^\t)) = K;(a;,FW(.,i,/j)), 

where, for every function G in L^(]R^"), the function Vci{.,G) is defined as in (1.1): 

(5.3) V,i{x,G) = V{x)+ j W{x-y)G{y,v)dydv. 
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One shows similarly that: 

(5.4) Vq{x,Ph{t)) = Vel{x,Uh{.,t)), 

where Uh is defined in (1.7). For all suitable functions A and B , let us denote by Mh{A,B) the Moyal 
bracket of A and B, defined in (3.2). With these notations, the function Uh{-,t) defined in (1.7) must verify: 

(5.5) Ou^^^^^^du^ ^ lM.(y.(n.(.,0),«.(.,0). 

For all functions A and B in C°°(IR^"), and for every integer k > 0, let Ck{A, B, .) the function defined in 
(3.4). We set Cq{A, B) = 0. One has Ci(A, B) = \{A, B}. We will choose the Uj in a such a way that the 
equation (5.5) is approximatively verified. 

The construction of the functions Uj of theorem 1.2 is detailed in the following proposition. 

Proposition 5.1. There exists a sequence of functions {X,t) — ?• Uj{X,t) on IR^" x IR (j >Q), such that: 

a) The function t Uj{.,t,h) is C°° on JR valued in W^°°'^(IR^"'), bounded in that space independently of h 
in (0, 1] and oft in every compact set of JR. 

b) One has: 

uo{X,i)) = FhiX), Uj{X,0,h) = 0, {j>l). 

c) For every N, the function un{X, t, h) verifies: 

(5-6) ^ + 2E^.^ = 7 E C,{v.i{;uA.,t,h)) , u,{.,t,h)). 

j=l ^ i+k+l=N+l 



In the sum (5.6) the indices j et £ are > and A: is > 1. 

Determination of uq. For N = 0, the equation (5.6) reduces to the Vlasov's equation: 

^+2E6-^ = E^^cK-o(.,o)^^ 

It is well known that there exists a unique solution of the Cauchy problem, verifying this equation and 
Uo{.,0,h) = Fh, where Fh is a given function in W°°^(5R^"'). One knows that then the function uq is 
continous on IR valued in H^°°^(1R^"). If Fh stays bounded in W°°^(]R^"') independently of h, it is also the 
case for uo{-, t, h). 

Determination of un (N > 1). For every N > 1, the equation (5.6) can be written as: 



Giv = ^ E Ck(Vci{uj{.,t,h)) , ue{.,t,h)y 



I 

j+k+e=N+i 
j<N,e<N 



One also requires that un^X, 0, h) — 0. To solve this equation, by dropping the parameter h for the sake of 
simplifying notations, let us denote by X — ^ <^t(-^) = {lit, X),p{t, X)) the the Hamiltonian flow, solution 
of: 

q'{t, X) = 2p{t, X) p'{t, X) = -VV{q{t, X)) - [ VW{q{t, X) - y)uo{y, r], t)dydrj, 
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verifying : 

q{Q,X) = x p(0,X)=e X={x,i). 

The function vn defined by VN{X,t) = UN{'^t{X),t) must verify: 



dvM dur, , , f 9W' 



1/ 



Gr,{X,t) = GN{MX),t). 

By using in the integral the change of variable {y,r]) = <pt{z,0^ whose jacobian equals 1, we see that 
must satisfy 

^{X,t) = GN{X,t) + [ A{X,Y,t)vN{Y,t)dY, 

eft 7]R2n 

Moreover, one must have vn{-,0) = 0. According to standard results on the Vlasov equation, one knows 



that Vuo(-, t) is in W (IR ), bounded when t varies in a compact set. The same is true for Vipt- If the Uj 
(0 < j < A'') have been built with the properties of proposition 5.1, one sees that GAr(.,i) is in H^°°^(]R^"), 
bounded when t varies in a compact set. It is also the case for GN{-,t). The resolution of the Cauchy 
problem verified hy vn and the one verified by uat is standard. 

□ 



To prove theorem 1.2, we will show that the functions F'^'^\t,h){X) defined in (1.11), starting from the 
Uj{.,t,h) of proposition 5.1, and the 
proposition is an intermediate step. 



Uj{.,t,h) of proposition 5.1, and the operators pl^\t) defined in (1.13), verify (1.12) and (1-14). The next 



Proposition 5.2. Let Ph{t) be an exact solution of (TDHF) verifying the hypotheses of theorems 1.1 and 
1.2. Let Hhit) he the operator defined in (4-3), where V{t) is the multiplication by Vq{x, ph{t)). Then the 
above constructed approximate solution p'j^\t) verify: 

(5-7) = [HH{t),P^^\t)]+OpZy\s\^\.,t)) 

where S\^\.^t) is in W°°^(S\?^) and verify, for every integer a: 

(5.8) ||V"5f )(.,OllLi(m-) < CMt)[h''^^ + h\\pH{t)-p^^\t)\\c^(n) 
where CaN{t) is a bounded function on every compact set. 

Proof. If the Uj are determined relying on proposition 5.1, the function f (-'^) defined in (1.11) verifies: 

(5.9) ^ + 2X:6^ = ^E'^'^'^(^«'(-'^^''H.,i,/^)) , F(^){.,t,h))+<^^''\.Xh), 

i=i fe=i 

where ^^^\.,t, h) is a function in W°°'^{M^'^), such that: 

(5.10) ||V«$W(.,t,/i)|Ui(]R2„) < h^'+^CMt)- 
We define an approximate version of the operator Hh{t) by setting : 

(5.11) ^r''(*) = -/i'A + Vd(FW(.,t,/i)). 
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Since F^^^ verifies (5.9), we may write: 



(5.12) 



ih- 



dt 



where the function Tj^^\.,t) is defined by: 

Tt''\.,t) = h^^^\.,t,h) + Ri''+'\v,,{.,F(^\.,t,h)),F^^\.,t,h)). 

(For all functions A et i? verifying the hypotheses of theorem 3.1, wc denote by R^^\A,B, .) the function 
associated by theorem 3.1. to such functions.) Then by the definition (5.3) of the map Vd, and point a) of 
proposition 5.1, we can write: 

(5.13) \\V''M■,F(''\.,t,h))\\L^^J^..) < C^Nit), ||V^^^W(.,t,/i))||ii(E2.) < CpN{t). 

Using these upper bounds and following theorem 3.1 on the Moyal bracket, we may write: 



V'Ri''+'\v4.,F(^H.,t,h)),F(^\.,t,h)) 



Li(]R2") 



Acccording to these upper bound estimates, and the estimates (5.10) of the derivatives of ^'^^\.,t,h), one 
has: 

(5.14) ||V"TW(.,i)||Li(iR-) < C\N{t)h''+\ 
According to (5.12) and since 

we write the equality (5.7) with: 

(5.15) 5f)(.,i) = TW(.,i) + M,(y,(.,p,(i))-F,(.,pi"''(t)),FW(.,t,ft)). 
One has: 



':{v,i.,p,{t))-V,i.,pi''\t))) 



Using all of these estimates and the norm estimates (5.13) of F'^^\.,t,h) and by using theorem 3.1 on 
the Moyal bracket, it results that: 



(5.16) \\w'^MH{vMt))-V,{pi''\t)),F(^\.,t,h))\\^^^^^^^<C^^^^^ 



The norm upper bound estimate 5.8) of S^\.,t) results from (5.15), (5.14) and (5.16). 

End of the proof of theorem, 1.2. Let Uh{t,s) and Gh{t,s) be the unitary propagator and the mapping, 
defined in proposition 4.1, associated to the operator Hh{t) of proposition 5.2. The comparison of equalities 
(4.21) (verified by the exact solution) and (5.7) (verified by the approximate solution ), and the Duhamel 
principle, allows us to write: 

(5.17) p,{t) - )(0 = Gu{t, s) {OpZ^\{S^^\., s)))) ds. 

Since Uh{t, s) is unitary, the map Gh{t, s) conserves the trace norm, and from this we may deduce that: 
l|pUi)-PrWll£MH) < IJ^ \\Opr\si''\.,s))\U^n) ds. 
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Using the upper bounds (5.8) of Proposition 5.2, and according to proposition 2.1, we obtain: 

\\Ph{t)-pi''\t)\\c^^'H) < ll^Cis)[h''+' + h\\p,{s)-pi''\s)\\c^^n)] ds. 
By the Gronwall lemma, we deduce, with a different constant that 

\\Ph{t)-pi''\t)\\c^n)<C{t)h''+\ 
Therefore point (1.14) of theorem 1.2 is proved. We deduce from this inequality and from (5.8) that: 

||V"4'^^(.,t)|Ui(R2„) < CMt)h''+^ 

where CaN{t) is a bounded function on every compact set. From proposition 2.1 and lemma 4.4, for every 
multi-index (a/3), the operators: 

/i-^-^ {ad Q{h)f{ad P{h)rGh{t,s){Opr''{si''\.,s)) 

are trace class, and their trace norm is independently boimded of (t, s) in a compact set of IR, and of h in 
(0,1]. By (5.17), for every multi-index {a, (3), and for every iV > 0, there exists a function CaN{t) > 0, 
bounded on every compact set of IR, such that: 

\\{ad Q{h)f{ad P{h)r[ph{t) - P^u\t))\\c^(H) < Caiv(t)/i^+^ 
Using proposition 2.1, point b), we deduce: 

(27r/i)-"||ar^'(p,(i)-pf)(i))|U.(^..)<C(t)/i(^+i)-^. 
In other words, with the notations of theorem 1.2: 

\\uh{.,t)-F^^\.,t,h))\\Li^^..) < C(t)h''. 
Theorem 1.2 results from the above considerations. 
Appendix A. Proof of proposition 2.1. 

The proof of proposition 2.1 calls upon a different notion of symbol. One can associate to every bounded 
operator A in ^ a Sh{A) sur IR=^" x M^" defined by: 

where the ^xh are defined in (2.14). An explicit computation of integrals shows that: 
(A2) I < > I =e-3tl^-^l' \\^xh\\ = l- 

Consequently: 



(^.3) Sh{A){X,Y) 



The interest of considering this function Sh{A) (which is up to a slight modification what FoUand [F] calls 
Wick symbol), is that it is related to the Weyl symbol, in both ways, by integral operators (proposition Al), 
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and that it can be majorized (in one norm) and minorized (in a different norm) by the trace norm of A 
(proposition A2). 

Proposition A.l. The Weyl symbol of an operator A is related to the function Sh{A) by 
(AA) Sh{A){X,Y) = (tt/i)-" / e-^(^-^)-(^-^) ar'''(^)(^) dZ, 

(^.5) cjZ''\A){Z) = 2"(27r/i)-2" / Sh{A){X,Y) Kh{X,Y,Z) dXdV, 

(A6) if^(X,r,Z) = e-*(^-^)(^-^)-^l^-^l'. 

By the definition formula (2.3) of the Weyl calculus, one has: 

A={nh)-^ [ i:zh a^''\A){Z) dZ, 

where T,zh is the operator defined in (2.4). An explicit computation shows that: 
(A 7) < ^zh'i'xh,'^Yh > ^ g_ 1 (z-x).(z-r) 

The equality (A. 4) follows. By the fundamental formula (2.16) of coherent states, on has: 
(^.8) A = (27r/i)-^" / < A<i/xh, > PxYhdXdY, 

where PxYh is the operator defined by: 

(^•9) PxYhf =< f, > "^Yh- 

One knows from (2.5) that: 

ar'^'iPxYhKZ) = 2"Tr{PxYh o Ez) = 2" < Ez^^-y^, ^xh > ■ 
By computation (A. 7) (where X et F are permuted ) we may deduce: 

ar"'(A)(Z) = 2"(27r/i)-2" / Sh{A){X, Y) \ < ^xh, ^Yh > \^e-i(^-^^^''-^UXdY. 

Using the equality (2.15) on the scalar product of coherent states, we obtain (A. 5). 
Proposition A. 2. Let A be a trace class operator and G a function in L^(]R^"). Then one has : 

(AlO) (27r/i)-2"^^^ I < A^xh,^Yh > G (^^) | dXdY < {2n)-^G\\L.^j^2nj\\A\\c^H), 



(All) WMc^-H) < (27r/i)-2" / 



<A^Xh,^Yh > 



dXdY. 
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Proof. Wo may write A = B1B2, where Bi and B2 arc Hilbert-Schmidt. By using the fundamental property 
(2.16) of coherent states one sees that, for all X et Y in J?^": 

< A^xh, "^Yh >=< B2^xh, Bl^Yh >= (27r/i)-" / uzhiX) vzh{Y) dZ, 

where we have denoted by uzh{X) =< B2^ xhi"^ zh > and Vzh{X) =< ^zh^Bi^^xh >■ Let Ih be the left 
hand side of (A. 10). By the Schur lemma: 

h < (27r/i)~^"/i"||G||j:,i(]R2n) / ||wz/i||i,2(]R2n) \\vzh\\L^{m?^)dZ. 
By (2.16), we have Wuzhh^M^n) = (27r/i)"/2|| 5**^^11 et \\vzh\\LHm.^'^) = (27r/i)"/2||Bi*zfc||. Hence : 

h < (27r/i)-2"/i"||G||ii(iR2„) / ||Bi*zft|| m^zhW dZ 
By the fundamental property of coherent states: 

(27r/i)-" / \\Bj^zhfdZ = {2nh)-'' [ < B* Bj^ zh,^ zh > dZ 

= TTiB*B,) = \\B,\\%^^^ 
where ||i3j||£2(-j^) is the Hilbert-Schmidt norm of Bj. Therefore: 

h < (27r)~"||G||i:,i(]R2„) ||Bi||£2(9^) ||S2||£2(-^). 

By taking the infimum over all the decompositions A = B1B2, one gets (A. 10). The inequality (A. 11) is 
then deduced from the equality (A. 8) since the operators PxYh have a trace norm equal to 1. 

□ 

Proof of proposition 2.1. For point a), we use the equality (A. 4), by integrating by parts , as in Rondeaux 
[R]. Hence we have shown: 

(27r/i)-2" / e-'^\''-^\' \Sh{A){X,Y)\dXdY <Ch-'' ^ /i^'"'+"''^/^||5^5|^^IUi(iR=-)- 

|a| + |/3|<2n+2 

One deduces then the upper bound estimate of the trace norm of A (point a) by using the second point of 
proposition A2. 

For points b) et c), we arc going to integrate by parts in the second equality (A. 5) of proposition A.l. One 
verifies that the function defined in (A. 6) is invariant by the following differential operator: 

L{h)Kh = Kh £(/i)= (^1+ '^"^'' ^ {1 + {Y-X)dx). 

Thus, equality (A. 5) implies, for every integer A'': 

\cjZ'\A){Z)\ < 2»(27r/i)-2" / \Kf,{X,Y,Z)\ \CL{h)f Sh{A){X,Y)\ dXdY. 
One verifies that: 

|i^,(x,y,z)| = e-*i^-^i^-^i^-^r 

19 



On chooses N = 2n + 2. There exists C > such that: 
\ar'iA)iZ)\ < ... 

...<C V h\"^/\2nh)-^^ [ \Kh{X,Y,Z)\G(^^^) \d^Sh{A){X,Y)\dXdY, 

|a|<2„+2 -^IR*" \ Vh J 

G(X) = (l + |X|)-2"-2. 

By the formula on the commutators: 
\ar'{A){Z)\ <... 

<C V h-^n-i\a\ + m/2 j ^-^\Z-^\--i,\X-Y\-Q(X-Y\ ^^^^^^^^^^^^^^^^^^^^ 
|a| + |/3|<2„+2 \ Vh J 

A^0h = iad P{h)T{adQ{h)fA. 
The above equahty can be also written as: 

\aZ'\A){Z)\ <... 

|a| + |/3|<2n+2 \ Vh J 

Inequality (2.13) is a consequence of the first point of proposition A. 2, and point b) follows easily. 
Appendix B. Proof of theorems 3.1 and 3.2. 

First step, common to both theorems 3.1 and 3.2. We know that, for all suitable functions F etG , Op^^^\F)o 
OpZ'\G) = OpZ'\C,{F,G,.), with: 

Ch{F, G, X) = (7r/i)-2" / e-^''^^-^'^-^^F{Y)G{Z)dYdZ 

where a is the symplectic form a{{x, (y, t])) = — xr). Consequently the Moyal bracket Mh{F, G, .) is 
defined by Mh{F,G,X) = Ch{F,G,X) — Ch{G,F,X). Thus, it suffices to write an asymptotic expansion 
Ch{F, G, .). We may write Ch{F, G, X) = 1), by setting, for every 6 e [0, 1]: 

e) = (nh)-^" I e-f<^(Y-x,z-x) p^Y)G{X + 0{Z - X))dYdZ. 
Consequently, for every integer N: 



JV-l 



Ch{F,G,X) = J2 l-,d!^^h{X,0) + Ri''\F,G,X) 



k\ 

k=0 



with: 



One sees, using integration by parts, that: 

k 



d^^X,0,h)= (^^^ {nh)-^'' j (,-^-K<Y-x,z-x) (a{Vi,V2)\F®G)yY,X + 0{Z-X))dYdZ. 
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If a function $ depends only on the Y variable one has (in the sense of distributions): 

and similarly if $ depends only on the Z variable. For ^ = 0, one has, by the above two equalities: 

d^^{X,0,h)=(^^^ a(Vi,V2)'=(F0G)(X,X). 

Therefore we do have indeed the equality (3.3) of theorem 3.1, by setting: 

{B.l) iiW(F,G,X) = ^W(F,G,X) - Ri''\G,F,X). 

It remains to obtain an upper bound for the norm of the two above terms. One also has: 

G, X) = (^) (7r/j)-2" / ii-3!l-lKf,{X, Y, ZMX, Y, Z, e)dYdZde 

where 

if;,(x,y,z) = e-T-(^-^.^-^) 

*(X, Y, Z, 0) = (a(Vi, V2)'^(i^ ® G)) (r, X + e{Z- X)). 
The function is invariant by the following operator: 

1 + + 4^^") (1 - h^y - hAz) 



Therefore, for all integers K and i: 

|V^4^)(F,G,X)| < (nh)-^- f (i_^l^|v^(*L)^*(X,y,Z,^)| dYdZdO. 
Consequently: 

{B.2) h'^'^\W^Rf\F,G,X)\<C M«+'3)/27„^(X,/i) 

a + P<l+2K+2N 
a>N,l3>N 

7„^(X,/i) = /i-2" / (l-ey"'-' ~^^ ] \V"F{Y)\\V^GiX+9{Z-X))\dYdZd9\ 

■/lR4"x[0,l] \ Vh J I 

End of the proof of theorem 3.1. We integrate the above equality (B.3) with respect to X, by making the 
change of variables: 

X = (1 - 6)-'^ {X-6Z) Y = Y Z = Z. 

We obtain : 

... < / (1 - ^)W-2n-i fi + \X-Y\ + \X-Z\ \ |v"i^(y)| |V''G(X)| dXdYdYdZdO. 

>/iR6»x[o,i] V vh J 
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If one has A'' > 2n + 1 and if one chooses K = 2n+l,we deduce, by using the Schur lemma, that: 

\\lM.,h)\\LiiM^.) < q|V"F|Up(iR2„) ||V^G|U,(]R2„). 
Adding these inequaUties, we obtain: 

o:>N,P>N 

By similarly proceeding for rI^\g, F, .), we arrive at the upper bound (3.5) of theorem 3.1. Point (3.6) is 
then deduced by proposition 2.1. 

End of the proof of theorem3.2. If one chooses = 2n + 1, it follows from (B.3) that 
By reporting in (B.2), then in (B.l), we obtain the majorization (3.7) of theorem 3.2. 

□ 
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